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C " !) ' Abstract. I review the status of research, conducted by a variety of independent 

^SJ ' groups, aimed at the eventual observation of Extreme Mass Ratio Inspirals 

(EMRIs) with gravitational wave detectors. EMRIs are binary systems in which 

. , one of the objects is much more massive than the other, and which are in a state 

3 , of dynamical evolution that is dominated by the effects of gravitational radiation. 

■^r • Although these systems are highly relativistic, with the smaller object moving 

' relative to the larger at nearly light-speed, they are well described by perturbative 

>D ' calculations which exploit the mass ratio as a natural small parameter. I review 

Cm the use of such approximations to generate waveforms needed by data analysis 

algorithms for observation. I also briefly review the status of developing the data 

CNj 1 analysis algorithms themselves. Although this article is almost entirely a review 

^ ' of previous work, it includes (as an appendix) a new analytical estimate for the 

1/^ ' time over which the influence of radiation on the binary itself is observationally 

negligible. 
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1. EMRIs and IMRIs 

This article is intended to be a brief review of a variety of works which share a 
common ultimate goal: the observation of Extreme Mass Ratio Inspirals (EMRIs) 
with gravitational wave detectors. I assume that the reader is familiar enough with 
the subject that she believes such work to be worthwhile, but not so familiar that she 
has kept up with the most recent progress (astrophysical motivation can be found in 
the introduction of almost any of the more recent references below). This article is 
intended to report on such progress. It is not intended to be an in-depth review, such 
as Refs. PI21|3]. With that said, I now very briefly describe EMRIs. 

An EMRI is made up of a non-spinning test mass /x orbiting a black hole with mass 
M and angular momentum of magnitude aM < M^ (I use units where G = c = 1). I 
will use Boyer-Lindquist coordinates (t, r, 9, (p) to describe the background geometry 
of the black hole. The test mass perturbs the spacetime metric gap so that it deviates 
from the Kerr metric g^'p'^{M, a) 

gaf3^g^'p'''{M,a) + ha0, (1.1) 

where the perturbation is of order the mass ratio h^p ^ fj,/M. Assuming that the 
test mass is restricted to a region near the boundary beyond which no bound Kerr 
geodesies exist (that is, near the radius nsco of the innermost stable circular orbit, 
T ^ nsco) the mass ratio fi/M determines the nature of the system's evolution. That 
is, it determines whether the test mass quickly plunges into the hole, or instead moves 
in some sort of persistent "orbital motion" . This can be seen from a crude scaling 
argument. Suppose that the test mass moves on some sort of bound orbit about the 
hole. A distant observer would measure the energy of the test mass to he E ^ fi and 
the period of the orbital motion to be Toib ~ M. The power of the observed radiation 
would hel dE/dt ^ (/x/Af)^. The timescale over which the orbital energy changes T^ad 
would be the ratio of the orbital energy to the radiative power Trad ~ NP / ji. The 
supposition that the test mass orbits the hole is then self-consistent when 

— ~ -^ « 1 . (1.2) 

31ad M 

The EMRIs with frequencies / ~ l/Toib in the band of LIS A- like detectors will 
have masses in the ranges 10^ < M/Mq < lO'^ and 1 < fJ,/MQ < 10^. To a good 
approximation, these systems satisfy the condition for orbit-like motion p.2|l . See 
Sec. Ill B of Ref. 4 for a more rigorous estimate of when to expect orbit-like motion. 
A close relative of EMRIs are asymmetric binaries with frequencies in the band 
of LIGO-like detectors. These intermediate mass ratio inspirals (IMRIs) have mass 
ranges 10^ < M/Mq < 10^ and 1 < /i/A/© < 10. IMRIs of course satisfy the orbit-like 
motion condition to a lesser extent than EMRIs, but it may still be useful to search for 
them with techniques that were designed for EMRIs. Searching for IMRIs with LIGO- 
like detectors has caught the attention of astronomers since it amounts to a search for 
intermediate mass black holes. It is attractive to those working with EMRIs because 
of the similar physics, and because searches can be performed immediately. LIGO's 
S5 science run could likely observe IMRIs out to a modest distance of roughly 10 to 
60 Mpc. Though this range is unlikely to yield detections, advanced LIGO will be 
sensitive to IMRIs out to a much more promising distance of about 0.2-0.9 Gpc 0- 

X This relation can of course be derived, but the sceptic might recall some other familiar example of 
waves where power is proportional to the square of the wave's amplitude. Here, the wave's amplitude 
is the metric perturbation h^yfj oc fi/M. 
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2. The EMR in EMRI 

It's useful to first become familiar with the "EMR" part of EMRIs before discussing 
the "I". That is, it's useful to understand the geodesic orbits of test particles bound 
to black holes before considering the effects of radiation. This is motivated in part 
because EMRIs will spend the majority of their lifetimes in the regime of orbit-like 
motion, where Eq. H1.2() holds. Also, recent developments in the description of these 
orbits have become powerful tools which should find applications beyond EMRIs. 

Compared to Newtonian orbits, strong field black hole orbits have unfamiliar 
properties. Orbits close to the hole have three distinct orbital frequencies. This 
is because, unlike weak-field orbits which are planar, the orbit is confined within a 
toroidal region with three degrees of freedom. One way to define the boundaries of 
this torus is to choose values for the three constants of geodesic motion: energy E, 
axial angular momentum L, and Carter constant Q (the Kerr-analog of the magnitude 
of the non-axial angular momentum). The following coordinate-based definition of the 
orbital torus is often more intuitive. 

As the orbit rotates azimuthaly about the spin axis of the hole, it bounces between 
two radii Tmin < r < r^^ax- The radial boundaries are often defined in terms of an 
eccentricity e and a semilatus rectum p, both of which conform to their Newtonian 
definition in the weak field 

^inin P ^inax P /"O 1 ^ 

~jr ~ i + e ' ~jr ~ i-e ■ ^ ' ' 

Some authors omit the mass M here, such that their p would have dimensions of 
length. The polar motion of the orbit will also be bounded by some minimum angle 
^min < 0- Since the black hole is symmetric under reflection about its equatorial plane, 
the other polar boundary is redundant 9 < t: — ^min- Alternatively, one can define 
these boundaries with an inclination angle§ t 

L + (sgn i)0,„i„ = I . (2.2) 

The sgn L term is in place so that i varies continuously from to 180° as orbits go 
from prograde to retrograde. For weak-field orbits, t is the angle between the orbital 
and equatorial planes. In the strong field however, orbits are not at all planar. In that 
case L is an indicator of the "polar thickness" of the orbital torus. It indicates only 
the boundary within which the orbit bounces in and out of the equatorial plane. 

Neither the azimuthal {(j>), polar {9), or the radial (r) motions are periodic 
functions of the time t measured by a distant observer's clock, or even of proper 
time T. The radial and polar motion are however periodic functions of Mino-time A 
(defined by dr/dX = r'^ + a'^ cos^ 9) 6 . If the orbit begins from r = r,nin and 9 = 6*111111, 
when A = 0, then 

00 00 

r(A) = ro+2 ^ r„ cos(nT^A) , 61(A) = |+2 ^ 9k cos(fcTeA) , (2.3) 

n=l fe=l 

where r„, 9k, and T^^e are constants. The functions t{X) and 0(A) have similar 
harmonic decompositions except that they each (i) increase linearly with A, and (ii) 
contain harmonics of both T^ and Tg. These decompositions have proved powerful 
tools for two reasons. First, as Fig.^demonstrates, the series H2.3|l generally converge 
rapidly. Even for orbits with large eccentricity or inclination, one need only compute 

§ This definition of t is different from the one found in Ref. 8 and elsewhere. In Ref. |Si, what I am 
currently calling l was defined as 6inc- 
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Figure 1. Fourier series coefficients for r(A) in the case of a black hole spin of 
a = 0.9M, and geodesies with p = 6, t = 40°. The analogous series for 6(A), (/'(A), 
and i(A) converge similarly. The two curves which level out (e = 0.7 and e = 0.9) 
appear to have exhausted the precision capabilities of the numerical method used 
to produce this plot. 



a small number of coefficients in the Fourier series in order to evaluate the orbit 
with great accuracy at arbitrary times. Second, and somewhat surprisingly, it turns 
out that by exploiting these decompositions, one can show that in the frequency 
domain associated with coordinate time t, most functions of these orbits have a discrete 
spectrum made up of frequencies [7] 



fmkn = mf^ + kfe + nfr 



(2.4) 



Here to, fc, and n are integers, while fr, fe, and /^ are fixed frequencies determined 
by the boundaries of the orbital torus. For example, if a test particle moves on a fixed 
geodesic, a distant observer will find that the resulting gravitational waves oscillate 
only at frequencies fmkn- 

In the presence of radiation, returning the "I" to EMRI, waveforms observed by 
distant observers would have a sliding-comb type of frequency-spectrum, a discrete 
spectrum which over time slides around on the frequency axis. Since these waveforms 
will likely be weak compared to detector noise, models for the evolution of these spectra 
will be needed in order to enable some form of matched filtering. The remainder of 
this article will review EMRI-related work which falls into one of two categories: (i) 
the calculation of waveforms, and (ii) the development of data analysis algorithms. 
Since my own research falls into the first category, I will spend more time on the 
discussion of waveform calculations than on data analysis development. Of course, 
my bias should not be taken as an indicator of relative importance. 
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3. Waveforms 

Efforts to compute EMRI waveforms can roughly be grouped into three categories: 
(i) Capra waveforms, (ii) Teukolsky waveforms, and (iii) kludge waveforms. I define 
the categories as follows (these definitions are described in the following subsections): 
If a waveform calculation is based on solving the MiSaTaQuWa equations |TO1 E] , or 
some higher order version of them, the result is a Capra waveform. If a waveform 
calculation is based on solving the Teukolsky equation [^ |^ |^, the result is a 
Teukolsky waveform. If a waveform calculation is based on a variety of formalisms, 
some of which may even make conflicting assumptions, the result is a kludge waveform. 
I have listed these categories roughly in order of increasing availability and decreasing 
accuracy. However, the ultimate categorization factor will be taken to be the method 
of calculation rather than either its accuracy or availability. I now describe these 
waveforms, and the status of efforts to compute them 

3.1. Capra waveforms 

The name Capra comes from the Hollywood film director Frank Capra (1897-1991) 
whose ranch (which now belongs to Caltech) served as the location for the first of 
a now-annual series of "Capra meetings" . The ultimate goal of these meetings is to 
compute a class of waveforms that incorporates the leading-order (in /i/M) effects 
of both the radiative and conservative parts of the gravitational self-force. The self- 
force acts on the test mass and is induced by that same object's perturbation of the 
background spacetime (it is a force in the sense that, instead of interpreting the motion 
of the test mass as a geodesic of some perturbed version of the background spacetime, 
it can be interpreted as being forced to deviate from geodesies of the background.). 
The radiative part of the self-force produces radiation at the black hole's horizon and 
at large radial distances. The conservative part does not produce such radiation, 
but nonetheless influences the world line of the small object, and the corresponding 
waveform |12| . Perhaps the crown jewels of efforts toward Capra waveforms are the 
equations of motion for a test particle under the influence of the flrst-order parts of the 
self-force in an arbitrary background spacetime. These are the so-called MiSaTaQuWa 
equations |ll)[ lll| . I categorize any waveform as a Capra waveform if its calculation 
was based on solving the MiSaTaQuWa equations, or some higher order version of 
them. Capra waveforms can be thought of as "holy grail waveforms". In terms 
of potential accuracy, they are the most ambitious. They are the only waveforms 
capable of being computed to an accuracy of order (/i/M)^. Realizing this potential 
will require computing all of the first-order parts of the self-force, and at least some 
of its second-order parts [OlITH] . 

The holy grail characterization of Capra waveforms applies both to their accuracy 
and availability. To date, no calculation of Capra waveforms exists (with the exception 
of analog problems for scalar or electromagnetic fields 52I)- Perhaps the most 
advanced effort toward Capra waveforms presented at the 2005 Capra meeting ^H] 
is the work by Barack and Lousto ^7]. They developed a numeric code that can 
compute a waveform for a test mass on an arbitrary worldline in a Schwarzschild 
background. The code evolves the metric perturbation in the Lorenz gauge, so that 
the input worldline can be given by the solution of the mode-sum representation of 
the MiSaTaQuWa equations j^ El 1201 (known only in the Lorenz gauge, but for 
any orbit in either Schwarzschild or Kerr backgrounds) . Instead of such a worldline. 
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a circular Schwarzschild geodesic is used in Ref . jlTj . They successfully compared the 
power radiated to infinity with results from calculations based on Teukolsky waveforms 
(discussed below) . Incorporating the solution of the mode-sum representation of the 
MiSaTaQuWa equations into the code from Ref. ^7] would produce the first of any 
subcategory of Capra waveforms. 

While progress toward Capra waveforms is steady, much work remains. Second- 
order versions of the MiSaTaQuWa equations and their solutions must be derived. 
Then, codes similar to the one used in Ref. ^^ must be developed before the full 
potential of Capra waveforms will be realized. More details on this subject can be 
found in Poisson's short manifesto [5l] or in his comprehensive treatise |5]. 

3.2. Teukolsky waveforms 

In 1972, Teukolsky derived a powerful equation for describing first-order radiative 
perturbations of black holes |22i i22j |24j . It applies to cases where the hole is perturbed 
by a scalar, neutrino, electromagnetic, or gravitational field. Although his is a partial 
differential equation, Teukolsky showed that it can be separated into a system of 
ordinary differential equations (DDEs). Solving the Teukolsky equation numerically 
is therefore far simpler than traditional numerical relativity 29 . I categorize any 
waveform as a Teukolsky waveform if its calculation was based on solving the Teukolsky 
equation. Since solutions of the Teukolsky equation describe only radiative first-order 
(in /i/M for the case of EMRIs) effects, Teukolsky waveforms have less potential 
accuracy than Capra waveforms which, at least in principle, can also describe second- 
order and conservative effects. 

The calculation of Teukolsky waveforms has been something of an industry since 
the early 1990's (see table I of Ref. "K"). These calculations exploit the condition 
that the test mass moves in an orbit-like fashion, so that Eq. (|1.2|l holds. The source 
term in the Teukolsky equation is taken to be a point-particle on a bound geodesic 
of the background spacetime. Since the equation's solution contains all the radiative 
information, it describes the waveform produced by the particle's motion, and also the 
effect of that radiation on the orbit. By iteration, the orbit evolves from a geodesic 
to an inspiral while the waveform evolves from only a "snapshot" of what an observer 
would see over a short time, to the entire EMRI waveform. This strategy has only 
recently been applied to generic black hole orbits with both eccentricity e 7^ and 
inclination t 7^ 0°, 180°. The radiation snapshots produced by these orbits have been 
computed numerically |8()l El , but the iteration needed to produce the full EMRI 
waveform is still under development. I now sketch a few of the details of these new 
developments. 

As is described in Sec. ^ many functions of generic black hole orbits have a 
discrete frequency-spectrum at frequencies /,„fc„, given by Eq. H2.4|l . The radiation 
produced by a test mass moving on that orbit is one such function. A distant observer 
would see that radiation as the following waveform 

h+-ih^=--Y,Y. E E ^^5z™fc„(0)e-'""'="(*-'-^+^'"^+*»"'=". (3.1) 

1=2 rn= — I k= — 00 n= — CO mkn 

Here h^ and hx are the two independent components of the metric perturbation, 
<^mkn = '^''^fmkm and both the Coefficients ^;^,j,„ and the functions 5';„ifc„(0) are 
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found by solving the ODEs that separate out from the Teukolsky equatfon ||. The 
phase constants Ximkn in Eq. H3.1|l are determined by the initial position of the test 
mass ^, and can be set to zero when considering a fixed orbit. Similar expressions for 
the radiative changes in the orbiting particle's energy E^ axial angular momentum L, 
and Carter constant Q have been derived PHI EH EH EZ| ■ For example, the average 
change in the orbital energy is given by 

-^)=2^ Z^ Z^ Z^ j:^ (^|^;mfen| + O^lrnkn \Zj^ikn\ J i (3-2) 

1—2 m— — l k—~oo n— — oo mkn 

where aimkn are simple constants that are analytically known |S], and Z^f.^ are 
similar to Zf^j^^. The geometry of the orbit dictates which terms are significant in 
Eqs. (|3.1|) and (|3.2|) . If the orbit is both circular and equatorial, all of the terms with 
/c 7^ and n ^ vanish. Similarly, for orbits with eccentricity above about 0.3 and 
inclination below about 60°, one captures better than half of the radiation by keeping 
only the terms with /c = |H]. 

Summing over only one of the terms in the parentheses of Eq. H3.2(l gives either 
the power radiated at large radial distances, or into the horizon. One normally would 
arrive at this formula by first deriving each sum independently, and then enforcing 
global energy conservation [SI]- This derivation has a major drawback. It is not 
generalizable to the evolution of the Carter constant. Alternatively, one can derive the 
same result as follows. First, solve the Teukolsky equation (this gives the radiation 
field, or equivalently Z^^f.^ and Z^f.^). Second, express the radiative self- force in 
terms of the radiation field. Third, express the average rate of change in the orbital 
constant of interest in terms of the radiative self-force. Galt'sov showed that, under 
general circumstances, this gives the traditional result [2] when applied to both energy 
and angular momentum |32j . At the time however, the validity of using only the 
radiative self-force was not known. Mino has since proved this to be valid 0, and 
the method has now been used to describe the evolution of the Carter constant in 
terms of the same quantities used for energy and angular momentum (namely Z^^.^ 
and 2'/^fe„) I13EH1ES1EZ1- Although the final equation for {dQ/dt) in Refs. [231^ 
appear different from the more concise result in Refs. |2HII27| . their equivalence has 
since been shown p^ . 

It has recently been demonstrated that, for an analogous case of an electric charge 
moving on a Newtonian orbit, one must account for the evolution of ximkn in Eq. H3.1|l 
in order to compute the full EMRI waveform to leading order in the charge to mass 
ratio P2I- Also, for the case considered in Ref. ^^, the evolution of ximkn was 
completely determined by the conservative self-force alone. Solutions to the Teukolsky 
equation contain only radiative information, so if this result carries over to the strong 
field regime relevant to EMRI observations, it will severely decrease the usefulness 

(I These quantities are defined in Sec. Ill of Ref. |H|. In the case of perturbations to a black hole's 
spacetime geometry, the leading order correction to the (otherwise vanishing) Weyl curvature scalar 
i/>4, a quantity that completely describes the radiation in the perturbed spacetime, satisfies the 
Teukolsky equation. In the case of perturbations from a bound test particle, 1/14 can be simply 
projected onto a basis of angular functions Simkni^) a-nd radial functions Rimknir) which both 
depend on the orbit's fundamental frequencies i^rnkn- The angular functions satisfy an ODE which, 
apart its dependence on tOrnkm is homogeneous. The radial functions satisfy an inhomogeneous ODE 
such that Rimknir ^ 00) = Zf^i^^f^{uj^kn,r) and Rimkn{r -^ r+) = Z,'^^^/°°(tj„fc„, r), where 
the event horizon is located at r = r-|- . The explicit (analytically known) definitions of the functions 
f^'°"{i^mknj r) are beyond the scope of this review. 
^ The exact dependence is shown in Eq. (8.29) of Ref. )28| . 
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of both Teukolsky waveforms and kludge waveforms (described below) . It should be 
emphasized however that the analysis in Ref. ^21 applies in the limit where gravity is 
Newtonian, and that some of their findings are expected to be atypical of the strong 
field regime (see the last several paragraphs of Sec. V in Ref. ^21)- Even if their results 
apply partially to the strong field, so that in the strong field the evolution of ximkn is 
observationally significant, but can be well approximated from the radiative self- force 
alone, the equations governing such an evolution have yet to be derived. Though 
such equations should in principle be obtainable from the analysis in Appendix C of 
Ref. ^1], the resolution to this problem remains a subject of current research |15|. 

There is also a history of efforts to solve the Teukolsky equation without fully 
exploiting its separability — only the ^-coordinate is separated such that the codes 
must solve a "2 -I- 1" dimensional partial differential equation. These "time-domain" 
calculations have yet to match the precision of their "frequency-domain" alternatives. 
In the most generic time-domain calculation to date (equatorial or circular Kerr 
orbits (31]) total fluxes of energy and angular momentum agreed with frequency- 
domain codes to within about 25%, whereas independent frequency-domain codes 
often agree to as many as six digits 0. However, recent improvements in time-domain 
methods, including the use of an adaptive mesh, are promising. Such codes have 
so far been applied to Schwarzschild jSS]- There agreement with frequency-domain 
codes improved by about a factor of ten over the previous generation of time-domain 
codes, agreeing to within about 0.01%. See Ref. ^5| for a discussion of the possible 
advantages of time-domain methods. 

3.3. Kludge waveforms 

1 categorize any waveform as a kludge waveform if its calculation was based on 
a collection of different formalisms, some of which may even make conflicting 
assumptions (i.e. solving the flat-space quadrupole formula for a particle on a 
relativistic black hole orbit). Kludge waveforms are in generally less accurate than 
Capra or Teukolsky waveforms. However as more rigorous waveforms have become 
available, kludge waveforms have made good on their promise of capturing the 
dominant features of the more realistic waveforms "Ad' . Kluge waveforms are readily 
available, and can be computed very quickly. This makes them the tool of choice 
when scoping out candidate data analysis techniques (discussed in the next section). 
Their flexible design allows one to experiment freely, although admittedly crudely, by 
adding and removing proposed physical effects which would be unimaginably difficult 
to incorporate into Capra or Teukolsky waveform calculations. For example, kludge 
waveforms are already available for speculative non-Kerr background space-times that 
may be used in straw-man/null-experiment tests of general relativity 37 38 . In 
principle, kludge waveforms can also include effects due to the conservative self-force 
|3()l I43| . Kludge waveforms may be lacking in rigor and perfection, but they overflow 
with availability and adaptability. 

The majority of kludge waveform calculations are numerical. I now describe the 
analytical exceptions. As early as 1963, Peters and Mathews derived waveforms for 
eccentric Newtonian orbits using the flat space quadrupole formula |39j . Barack and 
Cutler jini stiched together sequences of those waveforms by using post-Newtonian 
formulas to evolve the constants of orbital motion. This enabled them to account 
for relativistic effects including those due to the black hole's spin, precession of the 
perihelion, and Lense-Thirring precession. The post-Newtonian formalism is generally 
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Table 1 . Comparison of numerical solutions to the Teukolsky equation |S| 151 
with analytical post-Newtonian expansions of those solutions 1271 . Here the black 
hole has spin a = 0.9M, and the orbits have semilatus rectum p = 100. The 
dependence of the disagreement on the orbital parameters is the expected result 
of the hybrid calculation's use of a truncated Taylor expansion in both eccentricity 
e and inclination i. 



invalid for EMRIs, since the test mass moves at nearly light-speed. However, their 
waveforms were useful for making rough estimates of how well LISA could measure 
an EMRIs parameters 001 ■ These analytic kludge waveforms will also likely be used 
in the future for a LISA Mock Data Challenge, a project that will use simulated 
LISA data to test proposed data analysis techniques for a variety of sources (initial 
results will be announced at the 2006 GWDAW meeting). Another analytic effort, 
originally motivated as a tool for exploring convergence of post-Newtonian series 
and for improving waveforms for neutron star binaries, derives the post-Newtonian 
expansion to solutions of the Teukolsky equation . This hybrid technique assumes 
both a small mass ratio and slow motion. The most recent contribution to this effort 
applies to Kerr orbits which are both slightly inclined and slightly eccentric [27| , and it 
is anticipated that the same calculation will soon be complete for arbitrary inclination 
|44|. Table^shows that the hybrid calculation of radiative fluxes of energy and angular 
momentum j27) compared well with pure Teukolsky calculations [5j. 

The following is a rough outline of the procedure used to compute kluge waveforms 
numericaly: First, evolve the orbital constants by integrating equations of the form 



7C7 



(3.3) 



where £ = {E,L,Q), or f = {e,p,L), to obtain £{t). Since {E,L,Q) can be translated 
exactly to (e,p, i), and vice versa, the choice here is somewhat arbitrary. These 
equations might be, say, analytical post-Newtonian approximations, or some numerical 
formula gleamed from a fit to numerical data from solving the Teukolsky equation. 
Second, using these solutions, determine the inspiraling world line x(t) by integrating a 
system of geodesic equations of the form (e.g. Carter's equations 0J if the background 
is Kerr) 

^=V[£(t),a,M,/i]. (3.4) 

Third, obtain the waveform by solving an equation of the form 

ft-y(Xficld,t) = Wij(Xfield,x[i]) , (3.5) 

where here the subscript "field" denotes the position of the observer. Examples of such 
an equation might be the flat space quadruple, or quadruple plus octupole formulas. 
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Figure 2. Comparison of Teukolsky 'H"5' (black) and numeric kludge 1431 (red, 
or light gray) waveform snapshots for two different generic black hole orbits. The 
parameters of the black hole and the orbits are shown in the titles above the plots. 
This level of agreement continues for essentially all time. The short time range 
shown here was used so that the reader can get a feel for the overlap by eye. The 
overlap estimates are 99% for the top panel (over a time of 20,000 M) and 97% 
for the bottom panel (over a time of 15,000 M). 



This numerical kludge procedure has been followed by several groups using a 
variety of different approximate formulas at each step |42[ I36L US] . The first step of 
this procedure, integrating the equation that evolves orbital constants 1)3. 3|) . has been 
executed by Gair and Glampedakis for the case of generic black hole orbits [^ . In a 
separate calculation, for a fixed generic black hole orbit that does not inspiral, Babak 
et al 1481 obtained snapshots of generic kludge waveforms by completing steps two 
and three [solving Eqs. 1)3.4(1 and 1)3. 5|l with £ held fixed]. They found remarkable 
agreement when comparing with Teukolsky waveforms (see Fig. ^. Estimates of 
the overlap between those kludge waveforms 031 and Teukolsky waveforms O |5] 
were typically around 98%, for p > 6. This high level of agreement suggests that 
kludge waveforms may ultimately be used directly in EMRI detection algorithms, 
as opposed to being used only when exploring data analysis options (that is, data 
analysis analysis). These two works [201 E2] will likely soon be merged in order to 
produce generic kludge waveforms which represent the entire inspiral. 

Another recent set of interesting kludge waveforms are those produced by test 
particles moving on equatorial orbits of a "quasi-Kerr" background spacetime (a 
spacetime which is identical to Kerr, but which has a perturbed quadrupole moment, 
see Ref. ^j for further details). Those authors found that, by varying the orbital 
parameters, one could generally find a background-orbit pair (Kerr, geodesic orbit) 
which resulted in waveforms with a high overlap (> 90%) compared to some pair 
(quasi-Kerr, geodesic orbit). If this holds true in general, the common wisdom that 
the waveform is a unique signature of the black hole will be cast into doubt. It 
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seems likely that this "confusion problem" won't persist when one includes the effects 
of radiation and generic orbits. Radiation may shift the orbit-pairs away from the 
overlaping regime. Also, the high overlap required common orbital frequencies, but 
it is not yet clear that generic quasi-Kerr orbits have orbital frequencies (as opposed 
to continuous spectra). Nevertheless, the question of whether or not the confusion 
problem is general remains open. 

4. Data analysis 

Even if all three classes of waveforms (Capra, Teukolsky, and kludge) were readily 
available, the problem of extracting EMRI waveforms from data would remain non- 
trivial. In this section 1 briefly describe the status of efforts to develop EMRI data 
analysis algorithms for LISA. As promised, this section will be less detailed than the 
previous section on waveform calculations. A more thorough review of these issues 
can be found in Ref. 0S]. 

Gair et al jlH] estimate that a "brute force" matched filter search algorithm 
(akin to LIGO searches for neutron star binaries 021) that uses year-long template 
waveforms to find EMRls would require on the order of 10 templates, rendering 
it computationally impractical. Their favored practical alternative is a heirarchical 
search which begins with short duration templates (lasting a few weeks) of modest 
accuracy. It should be emphasized that the full year or so worth of data is used at this 
stage — it is only the waveforms which are short. One still needs a scheme for jumping 
from one short waveform to the next, but that scheme does not demand modeling 
the phase evolution for times longer than a few weeks. This segment of the search 
can be thought of as the "detection stage" in which EMRI candidates are identified, 
and relatively modest estimates of their parameters are made. The search would 
then proceed toward a "measurement stage" which focuses in on narrow regions of 
parameter space using increasingly accurate and longer-lasting templates, and which 
realizes EISA's full sensitivity (e.g. measuring /i, M, and a with fractional accuracy 
~ 10~* pn"). They estimate that such a scheme would yield anywhere from tens to 
thousands of EMRI observations over EISA's lifetime +. 

In the hierarchical search envisioned in Ref. |45j , Teukolsky waveforms (or perhaps 
even kludge waveforms) lasting up to a few weeks will likely suffice as templates for 
the initial detection stage, whereas Capra waveforms lasting up to a few years will 
be needed for the final measurement stage. Let's focus now on the detection stage. 
The number of waveform snapshots (waveforms produced by a single fixed geodesic) 
needed to produce sufficiently accurate detection waveforms is currently unknown. 
For example, it is possible that snapshots alone might suffice as detection templates. 
If so, then the tools for EMRI detection are available already |H]. The following 
simple scaling argument would seem to suggest that snapshots cannot be used as 
detection templates. Expand the <j) coordinate of the test mass' world fine (a quantity 
interchangeable with the phase of the waveform for purposes of this rough argument) 
as follows: 

<Pit) = (t>o + <j>Qt + -^ot^ + ■ ■ ■ . (4.1) 

Here (/)o, (j)o, and (/)o are constants. The first two terms would be exact for a circular- 
equatorial geodesic, while the third term represents the effect of radiation. A waveform 

+ This estimate neglects sources beyond a redshift 2 = 1. Although LISA could detect more distant 
EMRIs, population estimates at such distances are unreliable 1451 . 
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Figure 3. Snapshot timescale Tanap/M versus mass ratio fj,/M for a black hole 
spin of a = 0.5M, and a geodesic with e = 0.5, p = 10, and t = 0°. The (black) 
solid-line and circles are numerical estimates reported by Glampedakis and Babak 
1381 . The (blue) dashdot lines are the analytical estimate JA.IOI evaluated at the 
indicated orbital radii, r/M = 2, 4, 10, and 20. The (magenta) dotted line is the 
estimate from the scaling argument Tgnap = M^J Mj {i. If waveform snapshots are 
to be sufficient for detection, they must have Tsnap > 3 weeks. This boundary is 
shown as (red) dashed lines for a few different test masses /x. 



constructed from the geodesic terms alone would produce a phase error of order unity 



after a time Tgi 



where 
r2 



ioT, 



1. By dimensional analysis, one would expect 



■^ snap J ..iiv^iv^ vu^snap 

^qM ^ 1 and (fioM^ r^ fi/M. This suggests that the waveform snapshot would be 
valid for times shorter than 



Ts, 



.,„ap-A/VWM- (4.2) 

This equation predicts that, for EMRIs that could be observed with LISA, a single 
snapshot would accurately represent the waveform for times ranging from a few seconds 
to a day or so — hardly the few weeks demanded j32] of detection templates. However, 
this is only a scaling argument, and it neglects potentially significant coefficients which 
must diverge as the Newtonian limit is approached. Recently, Glampedakis and Babak 
estimated values for Tsnap by comparing snapshots of Teukolsky waveforms to kludge 
inspiral waveforms in the case of eccentric-equatorial orbits |38'. They found, even for 
orbits relatively close to the horizon, that Eg . I4.2l underestimates Tgnup by a significant 
factor Tsnap /(My^M/fi) ~ 100 (see Fig.|31l. An alternative analytical estimate (jA.lOp . 
shown in Fig. Eland derived in the appendix below, gives similar results. These results 
suggest that existing waveform snapshots will suffice for detection, at least for some 
subset of EMRIs that could be observed by LISA. 

Lastly, I want to address a topic which is relevant to EMRIs, but which is also of 
general interest for LISA data analysis as a whole. One of the key differences between 
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LISA and LIGO-like detectors is that LISA is anticipated to find an over abundance of 
sources. This very attractive feature also poses a potential problem in that waveforms 
from many different sources will overlap each other, complicating the task of cleanly 
identifying any one of the sources. For example, the waveforms produced by galactic 
compact binaries (GCBs, e.g. a pair of white dwarfs) are nearly monochromatic, so 
that a collection of thousands of these waveforms resembles a Fourier basis complete 
enough to fit essentially any smooth function in their frequency band. Therefore, 
an attempt to first subtract thousands of galactic binary waveforms from LISA data 
before searching for other signals is likely to (i) produce erroneous population statistics 
for the GCBs and (ii) doom further searches by washing away any remaining signals. 
In order to avoid this problem, LISA data analysis algorithms will simultaneously 
search for a variety of different sources. Algorithms that detect single waveforms in 
isolation are an essential ingredient to this process, however the isolated techniques 
will ultimately be merged into one "global fit" algorithm. 

When it comes to LISA's general need for global fit algorithms, EMRIs are no 
exception to the rule. Their lower frequencies will likely be hidden by a background 
of unresolved GCBs, and their observable band is expected to overlap with resolvable 
GCBs and mergers of massive black holes j321- Recent advances toward global fit data 
analysis techniques can be found in Refs. jlHlEni- Though these works do not deal with 
EMRIs, similar work must eventually do so, and these are good indicators of general 
progress toward global fit analysis for LISA. In Ref. 0H], a "reversible jump Markov 
chain Monte Carlo" technique is used to simulate the detection and measurement of 
100 monochromatic signals. In Ref. |4y| . similar techniques are implemented for up to 
10 GCBs. See also the two contributions by J. Crowder and by E. D. M. Wickham to 
this meeting (GWDAW 10). These works differ more in their numerical methods than 
in their theoretical foundations. They each represent implementations of Bayesian 
Inference |50j . 

5. Provocation 

The conference organizers asked for a review that would provoke discussion. In case I 
have failed to provoke anyone so far, this list of questions (and my guessed answers *) 
might: Could EMRI detections be made using only the existing waveform snapshots 
as detection templates? (This would work for some EMRIs with fi/M < 10~^.) Can 
kludge waveforms be used as detection templates? (Probably, but other waveforms will 
be needed in order to determine which kludges suffice.) Is the conservative self-force 
needed for detection? (No, but it will be needed for detailed EMRI measurements.) 
Are EMRI waveforms really unique signatures of the background spacetime, or is 
there instead a Kerr versus non-Kerr confusion problem? (Including the effects of 
radiation and generic orbits will show that there is no confusion problem.) What 
are the prospects for global fit data analysis techniques? (It is difficult to speculate 
quantitatively on the prospects for this field. However, this type of analysis is a hot 
topic in many fields that are otherwise unrelated to gravitational wave detection |50| . 
Due to the rapid growth in the field, as evidenced by the regular additions to lists of 
relevant papers [SJi it seems wise to develop as many independent global fit methods 
as possible.) Should we look for IMRIs with ground-based detectors? (Yes, little is 

* These guesses are for provocation purposes only. They should not be used in applications where 
injury or property damage may result if said guesses turn out to be wrong. 
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know about the sources, and searches can be done immediately.) Detailed analyses 
of any of these questions would likely make a significant contribution toward eventual 
EMRI observations. 
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Appendix A. Analytic estimate of Tgnap 

I thank Eanna Flanagan for permission to include this paraphrasing of his analytical 
estimate for the time Tgnap over which waveform snapshots are valid. 

Suppose that an EMRI's true orbital phase 0truc(i) can be approximated as a 
quadratic in time, as in Eq. H4.1|l . while the orbital phase used to compute a waveform 
snapshot template (/>tcmp(0 i^ only linear in time, as would be the case if the snapshot 
was produced by a test mass on a circular orbit. Take the true and approximate 
waveforms to have the form h = Acos[2(j){t)], where A is constant, and (p{t) is the 
corresponding orbital phase 

htvuc = AcOs(2<^truo) = ACOS(20O + 200^ + 4>ot'^) , (A.l) 

/itcmp = Acos(20tomp) = A COS (2 00 + 2(j)ot) . (A. 2) 

Also assume for simplicity that the average overlap integral over a time T is just a 
time integral, from — T/2 to T/2, of the product of two waveforms, divided by T. The 
average overlap of the waveforms /itruo and /itcmp is then 

^2 /.T/2 ^2 /.T/2 

Vavg{T) = T^ / dt COs(2(?!)truo + 20tcmp) + :^ / dt COs(20truo-20tomp) .(A. 3) 



^J J -T/2 ■^^ J~T/2 

The first term in this expression vanishes rapidly, and is presumably unobservable. 
The second term dominates since it approaches zero very slowly (typically on time 
scale that is many orders of magnitude longer than the time needed for the first term 
to vanish). The normalized overlap (in the sense that it goes to unity for a pair of 
identical waveforms) of these two waveforms is then given by 

1 /-^/^ 

Korm(T) ^- I dt COs(20truo " 2(/)tomp) ■ (A.4) 

^ J-T/2 
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Substitution of the phases read off from the waveforms (|A.1|I and HA.2|) gives 



Korm(T) = - / dt COs(0oi') = \'" , (A.5) 

where 

C{x) ^ I dt cos(7rtV2) , (A.6) 



^0 

is the Fresnel cosine integral, and where 




:,(r)=rW^. (A.7) 



As in Ref. |2H1) define the time at which the approximate waveform is invalid Tgnap to 
be the time at which the overlap Vnorm drops to 95%. Solving C{x)/x = 0.95 gives 
X ~ xq = 0.67 « 2/3, which gives 



TsU = y^cT • (A.8) 

Now use the Newtonian formula (from Sec. Ill of Ref. [S2| with J7 = 0o and (l — (J)q), 

4>o = ^-l^^My%'/' , (A.9) 

and approximate the orbital frequency 0o with Kepler's law r/Af = (Afc/jp)"^/'^, to 
find 



T 

-^ snap 



11/4 




(A.IO) 



As expected, this improved estimate diverges in the Newtonian limit r /M — > cxo. Over 
the range of innermost stable circular orbits of rotating black holes 1 < r/M < 9, it 
differs from the result of the simple scaling argument by a factor ranging roughly from 
38% to 160. 
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